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ABSTRACT 

I argue for two modes of gas giant planet formation and discuss the conditions 
under which each mode operates. Gas giant planets at disk radii r > 100 AU 
are likely to form in situ by disk instability, while core accretion plus gas capture 
remains the dominant formation mechanism for r < 100 AU. During the mass 
accretion phase, mass loading can push disks toward fragmentation conditions 
at large r. Massive, extended disks can fragment into clumps of a few to tens of 
Jupiter masses. This is confirmed by radiation hydrodynamics simulations. The 
two modes of gas giant formation should lead to a bimodal distribution of gas 
giant semi-major axes. Because core accretion is expected to be less efficient in 
low-metallicity systems, the ratio of gas giants at large r to planets at small r 
should increase with decreasing metallicity. 

Subject headings: planetary systems: formation — planetary systems: proto- 
planetary disks — hydrodynamics — instabilities — radiative transfer 

1. Introduction 

Both core accretion plus gas capture (e.g.. Pollack et al. 1996) and direct formation 
by disk instability (Cameron 1978; Boss 1997) can form gas giant planets in principle. In 
order for gravitational instabilities (GIs) to form gas giants directly, the Toomre (1964) 
Q = Csk/ttGE, where Cg is the sound speed, k is the epicyclic frequency, and E is the surface 
density, must approach unity and the local cooling time must be < the local orbit period 
(Gammie 2001; Rice et al. 2005; Rafikov 2005, 2007). An isothermal, \ow-Q is strongly 
susceptible to fragmentation (e.g., Tomley et al. 1994; Boss 1998; Nelson et al. 1998; Mayer 
et al. 2002; Pickett et al. 2003). 

Analytical work (e.g., Rafikov 2007) and radiation hydrodynamics simulations (Nelson 
et al. 2000; Cai et al. 2006, 2008; Boley et al. 2006, 2007b; Stamatellos & Whitworth 2008a; 
Boley & Durisen 2008) find that disk fragmentation inside r ~ tens of AU is unhkely because 
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regions with high coohng rates and low Q are rare (cf, most recently, Boss 2008 and Mayer et 
al. 2007). When Q is low, the gas is cool, so cooling times are typically long. When cooling 
times are short, the gas is hot, so Q is typically high. In contrast. Boss (2004) argues that 
convection can increase the cooling rates enough to cause fragmentation. Likewise, Mayer 
et al. (2007) report convection-induced fragmentation in their simulations when the mean 
molecular weight is increased from 2.4 to 2.lQ Boley et al. (2006) and Boley & Durisen 
(2008) argue in return that the reported convective motions are likely shock bores along the 
spiral shocks (Boley & Durisen 2006) or an artifact of suddenly changing the mean molecular 
weight. Moreover, Rafikov (2007) and Boley et al. (2006, 2007b) find that convection does 
not decrease cooling times enough to trigger fragmentation because sustained convection is 
regulated by the entropy gradient and energy must be radiated away at the disk's surface. 
Core accretion seems to be the only option for gas giant formation inside r ~ 100 AU. 
However, the recently-discovered substellar companions with wide semi-major axes a ~ 100 
AU and masses M ~ few to tens Mj (e.g., Luhman et al. 2006; Lafreniere et al. 2008; Kalas 
et al. 2008) challenge the interpretation that core accretion is the sole formation mechanism 
for gas giants. The susceptibility of disks to fragmentation at large r remains an open 
question. Stamatellos et al. (2007; hereafter SHW2007) and Stamatellos & Whitworth 
(2009; hereafter SW2009) find that extended, massive disks can fragment under realistic 
conditions (cf Boss 2006); but, their results are based on running simulations with highly 
unstable initial conditions, which as argued here, can significantly change the outcome of 
disk evolution. 

Consider a disk where the sound speed becomes constant with radius for large r. For 
S ~ r"'', Q ~ r'^Q ~ ^g-i s iQ^g as g < 1.5, the stabilizing shear contribution to Q falls 
off faster than the destabilizing self-gravity contribution, and the disk can become susceptible 
to GIs. Now consider the radiative cooling time for some fluid element trad = e/|V-F|, where 
e is the internal energy density of the gas. Using the free-streaming limit for the divergence 
of the flux, where | V ■ -F| = ApKpaT^ for low optical depth r = pupdz, the cooling time 
can be approximated by trad ~ Sc^/(4TcrT^), where p is the mass density at temperature 
T and Kp is the Planck mean opacity. For r <^ 1, radiative cooling becomes extremely 
inefficient. Because cooling in the optically thick disk is also inefficient, the r ~ 1 transition 
region may be the most conducive regime for disk fragmentation (e.g., Rafikov 2005). The 
cooling time in a disk with S ~ 1 g cm~^, r ~ 0.1, and T ~ 10 K, is trad ~ 600 yr, which is 
less than the orbital period around a 1 Mq star at r ~ 100 AU. 

What if the r ~ 1 and Q ~ 1 regions do not coincide? Mass loading, i.e., mass accretion 



^Boley et al. (2006) simulated a disk with a /x ~ 2.7; they did not see fragmentation. 
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onto the disk, is crucial for driving fragmentation at r > 100 AU. Numerical fragmentation 

studies have focused mainly on thermal energy balance arguments. However, there are 
additional considerations when evaluating the stability of a disk against fragmentation. The 
rate that Q changes in a local region of a disk is 

0? In Q 1 (i In (i In K, (i In E 

— = H . (1) 

dt 2 dt dt dt ^ ' 

When Gls are in a self-regulating phase, the effective Shakura & Sunyeaev (1973) a tends 
toward a few hundredths or less (e.g., Gammie 2001; Lodato & Rice 2004; Mejia et al. 2005; 
Boley et al. 2006), so the timescale for E and k, to change due to mass transport should be 
much larger than the local dynamical time. For fragmentation criteria studies of isolated 
disks, the thermal energy term is the most relevant. In order to avoid fragmentation, the 
local cooling time for the gas tcooi > /(t)^^, where is the local gas orbital period and 7(7) 
is some function of the adiabatic index of order unity. In the local model formalism, the ^cooi 
criterion is associated with a critical effective etc ~ 0.06 (Gammie 2001; Rice et al. 2005). 
This a-limit indicates the point at which local energy dissipation by Gls cannot balance 
cooling, and as a result, the disk fragments. 

In addition to energy transport, etc can be related to a critical mass flux Mc ~ STtacC^'E/D,. 
Substituting Q into this relation gives Mc ~ SacC^/ {GQ) (see also Goodman 2003; Matzner 
& Levin 2005). At large r, where the temperature of the disk approaches the envelope tem- 
perature, the ratio M^jM^ ~ 3a;c/<5, where Mg ~ c^/C is the mass accretion onto the disk 
from the envelope. For an envelope temperature of T = 30 K, M^. ~ 2 x 1O~^M0 yr~^. 
Because etc ~ 0.06, mass loading can operate at multiple times the maximum local transport 
flux. Episodes of non-local transport and/or disk fragmentation should occur for realistic 
envelope accretion rates. 

I present results from radiative hydrodynamics simulations that show mass loading can 
drive a disk toward fragmentation at large r. The fragments form clumps with masses greater 
than a few Mj. Throughout this Letter, I refer to the inner disk as the region inside r ~ 100 
AU and the extended disk as the region outside this boundary. This transition radius is 
illustrative, but represents a reasonable estimate for where fragmentation becomes possible 
(e.g, Matzner & Levin 2005). 



2. Methodology 

I use CHYMERA (Boley 2007) to model the formation of massive extended disks. 
CHYMERA is an Eulerian code that solves for the equations of hydrodynamics with self- 
gravity on an evenly spaced, cylindrical grid. The rotational states of molecular hydrogen 
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are taken into account (Boley et al. 2007a), and the radiative transfer scheme has passed a 
series of analytic tests. The radiative transfer algorithm and test results are presented in 
Boley et al. (2007b). D'Alessio et al. (2001) opacity tables are used for calculating optical 
depths, with a maximum grain radius amax = l/i^- In the following simulations, an equilib- 
rium ortho:parahydrogen ratio (0:P) is used. The gas is assumed to have a mixture of 0.73, 
0.25, and 0.02 for hydrogen, helium, and metals, respectively. 

Preliminary simulations using the CHYMERA radiative transfer scheme show that the 
radiative cooling time is shorter than the local orbital time and that the disk is effectively 
isothermal, in agreement with the estimate in §1. For modeling spiral waves, the normal 
radiative transfer scheme is stable in these simulations because the gas temperature is never 
far from the incident irradiation temperature. However, for some preliminary simulations, 
the cooling algorithm became numerically unstable immediately after clump formation due 
to the sudden increase in temperature, creating highly disparate radiative transfer and hy- 
drodynamics timescales. In order to study fragmentation at large r, the following radiative 
cooling algorithm is used for the simulations presented here. The divergence of the flux 



where A/V is the cell area-to- volume ratio, Tlrr is the incident irradiation on the disk at r, and 
fr = At + 1/Ar. The local optical depth across a cell is calculated by At — p(AtRosseiand(l ~ 
exp[-2Arpianck]) + Kpianck exp[-2Arpianck])V^/^, where Arpianck is an initial estimate using 
the Planck mean opacity. This cooling approximation goes to the free-streaming limit for 
small At and to zero for large At. In order to ensure that the algorithm is stable for large 
hydrodyamic time steps, the divergence of the flux is adjusted such that 



Here, e is the specific internal energy of the gas, Cequo is the internal energy of the gas if 
it were at Tjrr, Athydro is the Courant time step, and Atrad — pe/|V • F| is the radiative 
timescale. The irradiation temperature for these simulations is set to 30 K for all r. 

The central protostar's position is integrated sclf-consistently with a softened potential, 
where $star = GMstar/(|i" — r'p + s^)^^'^ for softening parameter s. In order to treat the force 
on the star, the mass in each cell is treated as a point mass at the cell's center, with the same 
softening parameter used for the star. The star's position is integrated from step i to i -|- 1 
by the following: = v^-^/^O.Sa^Af-^ v^+^Z^ = v' + 0.5a.'Af, and x^+i = x'-F v^+^/^Af . 
This algorithm is sufficient for maintaining, on average, the system's center of mass at the 
grid center. 



V-F^-{A/V)a{T^-T,t,)f;\ 



(2) 



y. ^adjusted ^ V • F exp[-(Aihydo/Airad)'] 

+ P(eequil - e)/Aihydro(l - exp[-(Aihydro/Airad)^])- 



(3) 
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Each disk is evolved on an r, (f) ,z — 256, 512, 64 grid, with a spatial resolution of Ar, 
rA(f), and Az = 2, 27rr/512, 2 AU. Mirror symmetry is assumed about the midplane, and the 
outer grid boundaries are outflow boundaries. There is an outflow boundary near the star, 
but negligible mass passes through it in these simulations. Mass is added to the grid near 
the top of the z outflow boundary between r = 60 and 300 AU. The added mass is given a 
density proflle of p ~ r~^, a specific angular momentum {G Mg^g^j-r^ / (r'^ + s^))^/^, an initial 
Vr = 0, and an initial Vz — — (2GMstar/^')^^^- Upon reaching a disk mass of 0.1 Mq, a random 
density perturbation is imposed with a maximum variation of ±10%. For all simulations, the 
disk is stable against GIs when the noise is added. For two of the simulations (see below) a 
softening of s = 20 AU was applied. This softening causes an error in the epicyclic frequency 
of just under 10% at r ~ 50 AU when compared with the Keplerian frequency and 3% at 
r ~ 100 AU. 



3. Simulations and Results 

Four simulations (SIMA, SIMB, SIMC, and SIMD) are shown in Figure 1. All disks 
except for SIMC fragment, and all disks evolve isothermally except in very high-density 
regions. Because the vertically integrated midplane r never becomes larger than a few, 
except in clumps, the above radiative transfer approximation is reasonable. Clumps reach 
temperatures in excess of 100 K. 

(1) SIMA: The protostar is set to 0.3 Mq. The initial mass loading ~ 1O~^M0 yr~^ 
until the disk mass — O.IMq, after which is reduced to 10~^Mq yr~^. The infall 
density profile is set to p = 1.5. A softening s = 20 AU is used because the accretion rate is 
less than what is used in the other simulations, requiring a longer evolution. The simulation 
is evolved for 16400 yr (3 P200, orbital periods at r ~ 200 AU), and reaches a disk mass of 
0.21 Mq before the simulation is stopped. Dense spiral waves develop, and a condensation 
forms near r ~ 70 AU. Fragmentation appears to be near wave corotation, in agreement 
with Durisen et al. (2008). By the end of the simulation, the clump has grown to 20 Mj 
and is located at r ~ 110 AU. As a cautionary check, a portion of this simulation was rerun 
with CHYMERA's normal radiative transfer algorithm. The disk behaved isothermally. 

(2) SIMB: The protostar is set to 1 Mq, and ~ IO^^Mq yr"^ for the duration of the 
simulation. The simulation is evolved for 5000 yr (~ 1.8 P200), and the disk grows to 0.52 
Mq. The infall density profile p = 1.5, and no softening is applied. Mass loading drives Q 
below unity outside r ~ 100 AU, and the disk fragments into 11 condensations by the time 
the simulation is stopped, ranging in mass from ~ 4 to 14 Mj. 
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(3) SIMC: Similar to SIMB, but the mass accretion is halted when = O.SSMq. At 
this mass, Q in a region near r ~ 90 AU drops below unity. This simulation investigates 
whether the disk can recover from the initial mass loading and avoid fragmentation. When 
Gls did set in, they rapidly transported mass both away and toward the star, so a softening 
of s = 20 AU was applied and the inner disk boundary was moved from r ~ 20 AU to r ~ 6 
AU. The initial simulation without softening shows a similar behavior, but was not evolved 
beyond the initial burst of Gls. The burst transports mass efficiently and Q reaches a mass- 
weighted average Q ~ 1-3 between r ~ 100 and 200 AU. This simulation is evolved for 9200 
yr (~ 3.3 P2oo)- Although a longer integration may be required to ascertain whether the disk 
will eventually fragment, the simulation does indicate that fragmentation is not guaranteed 
in the outer disk. 

(4) SIMD: The protostar is set to 1 M©, and the mass accretion rate is lO^^M© yr~^. 
As in SIMC, the mass accretion is halted once — O.33M0. The density profile for the 
infalling material p = 0.5. This choice for p places the minimum Q further out in the disk, 
such that Q first drops below unity near r ~ 140 AU. Unlike SIMC, the spiral waves are 
unable to redistribute the mass enough to avoid fragmentation. The disk forms a 6 Mj 
clump at r ~ 140 AU, which is subsequently transported to just inside r ~ 100 AU and 
grows to 11 Mj by the end of the simulation. The disk is evolved for 8600 yr (~ 3 P20o)- 

Wide semi-major axis, a, gas giants (WaCCs) and at least one brown dwarf (BD) are 
formed in these simulations. All disks go through strong burst-like phases, consistent with the 
arguments in §1. Although a high M is used for SIMB, SIMC, and SIMD, disk fragmentation 
is still expected for M ~ lO^^Mo yr^^, as seen in SIMA. The local Jeans length is resolved 
by at least four cells throughout the simulations (Truelove et al. 1997; Nelson 2006), where 
the local cell size is taken to be the geometric mean of the cell dimensions. These results 
indicate that fragmentation at large disk radii should be common and that WaGGs can be 
explained by in situ formation. These simulations demonstrate that the number or fragments 
and the fragment masses depend on how the disk is assembled. Studies investigating the 
clump mass spectrum of extended, massive protoplanetary disks, e.g., SW2009, are very 
sensitive to initial conditions. 

Although the fate of these clumps is unknown, the results of Mayer ct al. (2004) and 
SW2009 indicate that some condensations can survive. Radiative feedback, which is ne- 
glected in the current radiative cooling algorithm, should limit clump growth. Detailed 
simulations that can resolve the photosphere of a fragment are required to address this 
concern, which is beyond the resolution limit in these simulations. 

These disks are evolved through only one burst. Vorobyov & Basu (2005; hereafter 
VB2005) found in their 2D simulations that several bursts can occur during the disk accretion 
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phase. Whether a system forms WaGGs/BDs and retains them may be dependent on the 
details of the last GI burst. Because these simulations suggest that disks should be susceptible 
to fragmentation during their formation, they also support the possibility of a clump-driven 
outburst mechanism for FU Orionis objects (VB2005). The FU Orionis phenomenon is 
characterized by a rapid increase in the optical brightness of a young T Tauri object, typically 
5 magnitudes over a few to tens of years. Accretion from the inner disk onto the star 
is estimated to be as high as 10~^Mq yr~^ (Hartmann & Kenyon 1996). Because FU Ori 
objects have decay timescales of ~ 100 yr, an entire minimum mass solar nebula (~ 0.01 M©) 
can be accreted onto the protostar during an event. The best explanation for the optical 
outburst is a thermal instabihty (TI) (Bell & Lin 1994). However, prodigious mass flux 
through r ~ 0.1 AU seems to be required to drive the TI. A tidally-disrupted WaGG could 
supply the inner disk with 0.01 Mq mass, even when only the extended disk is gravitationally 
unstable. Helled et al. (2006) calculated the contraction time for a Jupiter-mass clump, and 
found that such a clump takes 3 x 10^ yr to reach a central temperature of 2000 K, i.e., 
the temperature required to dissociate enough H2 to lead to rapid collapse. This suggests 
that a WaGG has sufficient time to be transported inside r ~ 1 AU before it reaches 
mean densities that are high enough to avoid tidal disruption, where the disruption radius 
rt ~ 8 X 10~^(p/g cm-^)~^/^(M/Mo)^/^ AU. The Helled et al. model remains at a mean 
density of p ~ 2 x 10~^ g cm~'^ for ~ 10^ yr. This density corresponds to a disruption radius 
~ 1 AU for a 0.3 Mq star. An immediate observational signature of this mechanism would 
be a large shift in the radial velocity of the protostar just before an outburst. 



4. The Two Modes of Planet Formation 



Boss's (1997) advancement of the disk instabihty mechanism has spawned a decade's 
worth of work examining GIs as a formation mechanism for gas giant planets (see Durisen 
et al.'s 2007 review). One of the principal reasons for the mechanism remaining contested 
for so long is that, like core accretion, the mechanism works under the right conditions. 
Rafikov (2005, 2007), Boley et al. (2006, 2007b), Stamatellos & Whitworth 2008, and Boley 
& Durisen (2008) argue that radiative cooling timescales are too long for fragmentation out 
to r ~ 40 AU and that sustained convection does not cause fragmentation. This is also 
consistent with Nelson et al.'s (2000) 2D radiative hydrodynamics simulations, who assumed 
a polytropic vertical density structure in order to calculate a photosphere temperature. Gon- 
vection pushes the entropy gradient toward zero, so a vertical polytropic density stratification 
assumes efficient convection. Observationally, the planet-metallicity relationship (Valenti & 
Fischer 2005) indicates that planet formation favors a high-metallicity environment, which is 
strong evidence for core accretion (cf Boss 2005). The estimated core masses of Jupiter and 
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Saturn (Saumon & Guillot 2004; Militzer et al. 2008), along with the ice giants, support the 
core accretion mechanism. This evidence suggests that the dominant formation mechanism 
for gas giants inside r ~ 100 AU is core accretion. 

In contrast to the conditions inside r ~ 100 AU, optical depths should approach unity 
for a substantial Ar in the extended disk, as occurs in these simulations. Efficient radiative 
cooling, long orbital periods, and an equilibrium 0:P ratio combine to favor fragmentation 
in the extended disk. WaGG formation as a result of mass loading represents the first mode 
of planet formation, and takes place in the first 10^ yr of the disk's lifetime. Core accretion, 
which can continue after the main disk formation phase, represents the slower, second mode 
of gas giant formation. The survival of WaGGs will likely depend on the disk and mass 
accretion conditions during the last burst of GI activity. 

Both core accretion and disk instability are able to form gas giants under suitable 
conditions. In the inner disk, conditions favor core accretion, and in the extended disk, 
conditions favor disk instability. Scattering and planet-disk interactions should wash out 
any strict desert between the two formation regimes, but a bimodal population of gas giant 
planet semi-major axes should still be present. Because core accretion is expected to be less 
efficient in low-metallicity systems, the ratio of WaGGs to planets at small r should increase 
with decreasing metallicity. 

I would like to thank the referee for pointing out a complementary study by Rafikov 
(2009, astro-ph/0901.4739), which was posted during the review of this manuscript. Rafikov's 
work discusses fragmentation conditions and the stabihty of disks against high M's, while 
the study presented here demonstrates that a disk can be pushed realistically toward frag- 
mentation. 

I thank R. H. Durisen, L. Mayer, G. Lake, R. Teyssier, O. Agertz, D. Stamatellos, and 
the referee for comments that improved this manuscript. This research was supported by a 
Swiss Federal Grant and the University of Zurich. The presented simulations were run on 
NASA Advanced Supercomputing facilities. 



-9- 

REFERENCES 

Bell, K. R., & Lin, D. N. C. 1994, ApJ, 427, 987 
Boley, A. C. 2007, Ph.D. Thesis, Indiana University 
Boley, A. C, & Durisen, R. H. 2006, ApJ, 641, 534 
-. 2008, ApJ, 685, 1193 

Boley, A. C, Durisen, R. H., Nordlund, A, & Lord, J. 2007b, ApJ, 665, 1254 

Boley A. C, Hartquist, T. W., Durisen, R. H., & Michael, S. 2007a, ApJ, 656, L89 

Boley, A. C, Mejia, A. C, Durisen, R. H., Cai, K., Pickett, M. K., & D'Alessio, P. 2006, 
651, 517 

Boss, A. P. 1997, Sci, 276, 5320, 1836 

-. 1998, ApJ, 503, 923 

-. 2004, ApJ, 610, 456 

-. 2005, ApJ, 629, 535 

-. 2006, ApJ, 637, L137 

-. 2008, ApJ, 677, 607 

Cai, K., Durisen, R. H., Michael, S., Boley, A. C, Mejia, A. C, Pickett, M. K., & D'Alessio, 
P. 2006, ApJ, 636, L149 

Cameron, A. G. .W. 1978, M&P, 18, 5 

D'Alessio, P., Calvet, N., & Hartmann, L. 2001, ApJ, 553, 321 

Durisen, R. H., Boss, A., Mayer, L., Nelson, A., Quinn, T., & Rice, K. 2007, in Protostars 
and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil (Tucson: Univ. Arizona Press), 
607 

Durisen, R. H., Hartquist, T. W., & Pickett, M. K. 2008, Ap&SS, 317, 3D 

-. 2001, ApJ, 553, 174 

Goodman, J. 2003, MNRAS, 339, 937 



Hartmann, L., & Kenyon, S. J. 1996, ARA&A, 34, 207 

Helled, R., Podolak, M., & Kovetz, A. 2006, Icarus, 185, Issue 1, 64 

Kalas, P., Graham, J. R., Chiang, E., Fitzgerald, M. P., Clampin, M., Kite, E. S., Stapelfeldt, 
K., Marois, C, Krist, J. 2008, arXiv:0 811.1994| 

Lafreniere, D., Jayawardhana, R., & van Kerkwijk, M. H. (2008), ApJ, 689, L153 

Lodato, G., & Rice, W. K. M. 2004, MNRAS, 351, 630 

Luhman, K. L., Wilson, J. C., Brandner, W., Skrutskie, M. P., Nelson, M. J., Smith, J. D., 
Peterson, D. E., Gushing, M. G., & Young, E. 2006, ApJL, 649, 894 

Matzner, G. D., & Levin, Y. 2005, ApJ, 628, 817 

Mayer, L., Lufkin, G., Quinn, T., & Wadsley, J. 2007, ApJ, 661, 77 

Mayer, L., Luflkin, G., Quinn, T., & Wadsley, J. 2002, Sci, 298, 1756 

Mayer, L., Quinn, T., Wadsley, J., & Stadel, J. 2004, ApJ, 609, 1045 

Militzer, B., Hubbard, W. B., Vorberger, J., Tamblyn, I., & Bonev, S. A. 2008, ApJ, 688, 
L45 

Mejfa, A. G., Durisen, R. H., Pickett, M. K., & Gai, K. 2005, ApJ, 619, 1098 
Nelson, A. F. 2006, MNRAS, 373, 1039 

Nelson, A. P., Benz, W., Adams, P. G., & Arnett, D. 1998, ApJ, 502, 342 
Nelson, A. P., Benz, W., & Ruzmaikina, T. V. 2000, ApJ, 529, 357 

Pickett, B. K., Mejia, A. G., Durisen, R. H., Gassen, P. M., Berry, D. K., & Link, R. P. 2003^ 
ApJ, 590, 1060 

Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, J. J., Podolak, M., & Greenzweig, 
Y. 1996, Icarus, 124, 63 

Rafikov, R. R. 2005, ApJ, 621, L69 

2007, ApJ, 662, 642 

Rice, W. K. M, Lodato, G., & Armitage, P. J. 2005, MNRAS, 364, L56 

Saumon, D., & Guillot, T. 2004, ApJ, 609, 1170 



Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Stamatellos, D., Rubber, D. A., & Whitworth, A. P. 2007, MNRAS, 382, L30 
Stamatellos, D., & Whitworth, A. R 2008, A&A, 480, 879 
Stamatellos, D., & Whitworth, A. P. 2009, MNRAS, 392, 413 
Tomley, L., Steiman-Cameron, T. Y., Cassen, P. 1994, ApJ, 422, 850 
Toomre, A. 1964, ApJ, 139, 1217 

Truelove, J. K., Klein, R. 1., McKee, C. F., HoUiman, J. H., II, Howell, L. H., Greenough, 
J. A. 1997, ApJ, 489, L179 

Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141 

Vorobyov, E. I., & Basu, S. 2005, ApJ, 633, 137 



This preprint was prepared with the AAS IATeX macros v5.2. 



- 12 - 




Fig. 1. — Surface density snapshots at the end of each simulation. Strong bursts of GI activity 
occur in each disk. For three of the simulations, these bursts lead to fragmentation and clump 
formation. Movies of these simulations can be viewed at |http: / /www.aaroncboley.net, under 
the Movies tab. 



